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Transverse momentum fluctuations can be understood as resulting from clustering of strings or 
partons. Data allows to distinguish clustering without percolation, from clustering with percolation. 
Percolation is clearly favored by data. 
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Event by event fluctuations of total transverse momen- 
tum, P, in high energy hadronic and heavy ion coUisions, 
may give relevant information on the thermodynamical 
properties of dense hadronic matter, and, in particular, 
on the existence of a phase transition P, . The prob- 
lem has been discussed in a series of experimental papers 

In this paper we argue that the occurrence of P fluc- 
tuations is not necessarily a signal of a thermal phase 
transition. It may rather be a geometrical effect, due to 
clustering and percolation of strings 0, 0] . 
Let us suppose that in hadron-hadron and nucleus- 
nucleus collisions one always produces a cluster of N in- 
dependent identical strings. For the P distribution one 
can then write 



N 

f(p)N= I ua{p,)idp, s{p-Y,p^) 

i—l 



From Eq. ^ it follows that 

{P)n^N{P)i 

and 



{P^)n=N{P^),+N{N-1){P)1 



such that 



{P')^-{P)%=N[{P^),-{P)l)] 



(1) 



(2) 



(3) 



(4) 



as expected for independent strings. 

However, if one considers fluctuations in the number N 

of strings per cluster, one obtains 

{P')-{P)' = [{N') - (N)^] {P')i + {N) [{P''), {P)l] 

(5) 

Note that if one assumes that particles from the single 
string are emitted independently (Poisson like distribu- 
tion), one arrives at: 



{P\-{P)l=n{{p')^{pr) 



(6) 



where n is the string average multiplicity and p refers to 
the particle transverse momentum. 

Introducing now the quantity Fp^ (see, for instance, 
Refs. 3, 4]), 



(7) 



with = {P-{P)f, = (j)-{p)f and (n) = {N)n, 
such that Fp^ = for random emission, we obtain 



where 



1 _ (iV2) - (Ar)2 
k ^ 



(8) 



(9) 



is the normalized variance of the cluster distribution, and 



(10) 



is characteristic of the string particle production mech- 
anism. The inclusion of color summation factors when 
forming clusters of strings, not considered here, may af- 
fect the value of c, Eq. H1U|I . In any case, physics is 
mostly determined by {N)/k. 

Experimentally, data on fp^,, Ref. Q, show that this 
quantity increases as matter density (number of partic- 
ipant nucleons) increases up to some point starting to 
decreases for larger densities. Making use of two very 
simple string clustering models. Model I without perco- 
lation and Model II with percolation, we show that only 
Model II, with percolation, is consistent with data. 
Let us consider the (two dimensional) problem of dis- 
tributing Ns strings with a definite transverse area, in a 
total area M (the impact parameter transverse interac- 
tion area). A relevant parameter is the transverse area 
density rj: 



*Also at Escola Secundaria da Ramada, Portugal. 



M' 



(11) 
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Let us now assume that in a collision clusters of strings 
are formed, P{N) being the probability of having a iV- 
clustcr(cluster with N strings). 

In principle - by using Monte Carlo simulations - it is 
possible to calculate {Nc), the average number of clusters, 
(N) the average number of strings per cluster, higher 
moments of P{N), and the average area occupied by the 
clusters (A). 

There are two sum-rules which have, in general, to be 
satisfied: 

i. {N,){N)^N,, (12) 
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{N,){A)=M{l~e-^). 



(13) 



The first sum-rule expresses the conservation of number 
of strings and the second is simply the area sum-rule: 
the right hand side of this sum-rule is the overall area 
occupied by the clusters, Refs. Note that Eq. l(T^ 

and Eq. H13|) are invariants, independent of percolation. 
Percolation affects (A), (N) and {N)c, but not Ng and 
M{l-e-'>). 

In the next two sections we introduce the two Models 
that illustrate our arguments. 



Model I - No percolation 

This model is essentially equivalent to the problem of 
distributing Ng coins among M boxes. The total area is 
M, each box has area 1, and (A) — 1. From Eq. l(T^ 
and Eq. (|13|l one immediately obtains 



{N) = 



(14) 



The second moment (A^^) of P{N) can also be computed 
and we obtain for k the following equation. 



(iV2) _ (7V)2 l-(l + 77)e-'' 



(15) 



Note that as 77 — + 0, (N) — > 1 and k ~> 00. This means 
that when the coin density is very low, the clusters have 
just one coin and the distribution is a 5-function ((iV^) — 
(TV) 2 ^ 0). In the ?7 ^ 00 limit, {N'^)-{Nf -> (N) - ry, 
increases with 77. The behavior of the quantity of interest 
to us, {N)/k, is show, in the case of Model I, in Fig.(l). 



Model II - Percolation 

In the present context percolation means that the av- 
erage area of the cluster can not be constant: it must 
increase with density. In the low density limit, we can 
accept that Model I is correct, i.e.. 



(A) 



1. 



(16) 



FIG. 1: Average number of strings per cluster multiple by 
the normalize variance of the cluster distribution as function 
of dimension less density 77 in the context of AA collisions 
geometry, for Model I and Model 11. In both Models M = 16. 
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FIG. 2: Avarage area occupied by a cluster as function of 
dimension less density rj for AA collisions geometry, for Model 
I and Model II. In both Models M = 16. 



However, in the rj 00 limit, percolation requires 



M(l-e-"), 



(17) 



or, (Nc) I. In this case, P{N) is dominated by big size 
clusters. Percolation implies that we have a 5-function 
distribution, not only in the 77 — limit, P{N) — 6{N — 
1), but a ^-function distribution, as well, in the 77 — > cx3 
limit, P{N) -> d{N- (N)). 

In order to satisfy Eq. (|16|l and Eq. H17|l we choose for 
(A), as a function of 77, the function 



(A) = /(77) [M(l 



-1]+1 



(18) 



where 7(77) represents a percolation function, 7(77) = 
means no percolation, (A) = 1, f{r]) = 1 means full 
percolation, with the Heavyside function approximated 
by a smoother function, because of the finite size of the 
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systems in simulation, (see last paper in Ref. (g) 



fiv) 



(19) 



with this function we fit the fraction of occupied area 
by all string in the transverse plane and we obtain for 
a ~ 0.85 ± 0.01 and rjc ~ 1.155 ± 0.015 this parameter is 
the percolation threshold. 

In Fig. (2) we show (A) as a function, of ry for Model I 
and II. From Eqs. (fT^ . (fT!^ and ifTI^ one immediately 
obtains 

{N) = [Mil - e-") - 1] + 1) (20) 

with limits (N) —> 1 as 77 ^ 0, as in Eq. p4|l . and 
(N) Ns, as ?] 00 (the average cluster, asymptoti- 
cally, contains all the strings). 

In order to estimate the second moment (iV^), or the 
variance, one has to keep in mind that: 

(N^) - (N)^ > , 



(N') {N) 



2 '1^0 



0, 



(21) 
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By making use of Eqs. (|21() we write 



-(l+»7)(l-cxp(-,7))] 



1 - e-^ 



(22) 



It is easily seen that Eq. (|22|) satisfies the constraints of 
Eqs. PT|). 

In Fig.(l) we also show the quantity {N)/k for Model II. 
It is clear what is the essential difference between Model 
I (No Percolation) and Model II (Percolation). While in 
Model I, {N'^) — {N)'^ always increases with ij, thus origi- 
nating increasing P fluctuations, in Model II, as 77 ^ cxd, 
{N'^) — {Ny starts to decrease, thus originating decreas- 
ing P fluctuations. 



Data compared with models 

Before comparing Models I and II with data, by mak- 
ing use of Eq. ^ , we would like to give estimates for the 
values of M and c. 

In string percolation the variable 77, Eq. ((TT|l . is defined 

as 



(23) 



where r is the transverse radius of the string (we assume 
r = 0.2 fm) and R is the average interaction radius. It is 



^ 
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FIG. 3: Data from PHENIX experiment compared with 
Model I(M = 220, c = 0.220, = 8) and Model II(M = 
420, c = 0.430, iVf = 8) as a function of Npart- 



natural to identify Af , which has the meaning of an area, 
with [R/r)"^. We then have 



M < 



Ra 
r 



(24) 



where Ra is the radius of the nucleus with, for gold, 
Ra = 6.6 fm. We thus obtain the bounds 



25 < M < 1100, 



(25) 



where the left hand side represents the limit of pp colli- 
sion. Concerning the parameter c, Eq. (|10|l . the quan- 
tity {pf/iip'^) - (pf) = 1/(1 - 7r/4), if we admit, as in 
Schwinger model, that the particle px distribution of the 
string is gaussian. The quantity n, the string particle 
density, was in Ref. ^3 estimated to be 0.7, for valence 
strings, and 20 times smaller for sea strings. We obtain 
the bounds. 



0.163 < c < 3.262 



(26) 



Note that c should be a decreasing function of energy, as 
the fraction of sea strings increases with energy. 
The number of produced strings is obtained from the 
number of participants by making use of the approximate 
usual formula 



N, = NfN 



4/3 



(27) 



where Nf is the average number of strings in pp collision, 
NP is assumed to be equal 8 at y/s — 200 GeV, which 
agree with Ref. and Na — Npart f^- 
In Fig. (3) we compare our Models I and II with data 
published in Ref. The parameters M and c were, 

within bounds, inequalities Eq. I|25|) and Eq. (|26f) . ad- 
justed separatly for each one of the models. 
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In conclusion, data clearly favour Model II, with perco- and Roberto Ugoccioni for advice and discussions, 
lation. 
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